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The theory of extremely low frequency radio wave propagation from vertical ami 
horizontal electric dipoles in a half space, separated by an infinite slab from another half 
space, is discussed and application is made to the specific case of the sea water-atmosphere- 
ionosphere problem, with dipoles located in the sea water. Each of the media is assumed 
homogeneous and isotropic. When attention is restricted to the frequency range 1 to 1,000 
c/s, integration in the complex plane leads to consideration of the pole corresponding to the 
TEM mode of transmission and two branch cut integrals. One of these (that giving rise to 
propagation of energy along and in the ionosphere) is found to be important in tlie case of the 
horizontal dipole. 

1. Introduction 

Tlie ^^mode tlieory'^ of VLF and ELF radio wave propagation lias been well developed by 
several investigators (such as for example Wait [1960b] and Brekhovskikli [1960] who also give 
extensive bibliographies). An integral representing the Hertz potential in any laA^er of a 
multilayered mediiini due to a source in the same or an}^ other of the layers, is readily obtainable, 
the most common representation being as a superposition of cylindrical wave functions. 
Typical integrals of this sort are not easily evaluated, however, except in the very simplest cases. 
For this reason, most of the analytical effort in tjie past has been confined to models associated with 
rather common physical situations, as for example the propagation of seismic waves in the 
earth [Pekeris, 1948], or the electromagnetic wave radiation from lightning strokes in the 
atmosphere [Schumann, 1952]. 

In this paper attention will be directed to the particular three-layer model which simulates, 
in plane parallel geometry, the sea water-atmosphere-ionosphere problem of radiation from 
electric dipoles submerged in the sea water. The situation is illustrated in figure 1. It is 
expected that the results will approximate those of a spherical model reasonably well out to 
distances from the antenna on the order of the earth ^s radius [Wait and Carter, I960]. The 
equivalent two layer problem (dipoles submerged in sea water with the atmosphere assumed 
homogeneous and semi-infinite) was first investigated by R. K. Moore [1951] and later by 
Bafios and Wesley [1953, 1954]. 

As is customary in such problems, integrals expresssing the Hertz potentials are deformed 
into the complex plane and their determination is reduced to the evaluation of residues at 
poles, plus integrations around appropriate branch cuts. In the frequency range 1 to 1,000 
c/s only one of the poles — that corresponding to the TEM mode of transmission— is of impor- 
tance at appreciable distances. Asymptotic estimates of the branch line integrals are made, 
and it is found that the contribution from one of these, in the case of the horizontal dipole 
and the lower portion of the frequency range, is not generally negligible in comparison with 
the mode solution. This corresponds physically to the fact that appreciable energy is prop- 
agated along and within the ionosphere. 



• This work was performed at the University of New Mexico under Nonr contract 2798(01). 
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Figure 1. Geometry of the problem. 



2. Development of the Hertz Potentials 
2.1. Vertical Electric Dipole 

The development in this and succeeding sections is given in terms of the Hertz potential, 
designated by the vector tt. The boundary conditions can be satisified, for the vertical dipole, 
with a single component, which will be considered a ^-component, corresponding to dipole 
orientation along the s-axis of a cylindrical coordinate system. Subscripts 1, 2, and 3 will 
later be identified with sea water, atmosphere, and ionosphere in that order. Each of the 
three media is assumed homogeneous and isotropic. The time dependence is taken to be 
as e'''^ and is thereafter suppressed. 

The parameter kn is defined by 

kl = CD^IJL„e^—ioOfln(^n, 71=^1,2,3 

with the understanding that the square root having negative imaginary part is to be taken. 
This may also be written 



kl=o)'^jd^'eri, 



71=1,2,3 



where €n=en+^ 



is the so-called ''complex permittivity." The quantities €,„ fi^, and a^ are permittivity, per- 
meabihty, and conductivity, respectively. M. K. S. units are used throughout. 

The water-air boundary is designated by the plane z=0, and the air-ionosphere boundary 
by z=h. The dipole location is at z= — 2q (see fig. 1). Representations for the Hertz po- 
tentials are made as superpositions of cylindrical wave functions, and upon application of the 
boundary conditions we find for regions 1 and 2 (where we are mainly concerned with the fields.) 
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Ks+K2 tanli (7:/32A) 



K2+KZ tanh (^M) 



in whichi^„(X)=--^,n=l,2,3 

witli /?;, a separation parameter defined by 

and Pn=(Jn+i^en' The inia^'inary part of P^ is taken as negative. 
The identity 

j,{\p)=\[Tn'\\p)+Hi?\\p)] 

allows the integrals for ttj and ^2 to be written as the snm of two integrals of which the first is 
deformed into tlie first qnadrant of the complex X-plane; the second into the fourth quadrant, 
as shown in figure 2. It can be shown that the only contributions arise from integrations 
around branch cuts 2 and 3, and the resichies at poles. At appreciable distances the contribu- 
tion from 1)raiicli cut W will be small, so that it can ])e omitted from consideration. 



Figure 2. Contour integration in the 4^h quadrant of the 
complex \-plane. 
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The formal expression for the residues at the poles is, for region 1, 

-^ 2: ^(,, ^^^^^ ^^ g ■ ; (2.3) 

where Pi^^{k\—\f)^''^, and the Xj are the zeros of (i^i+>2'2). 
For region 2 the residues at poles are 

2.MJZ H^^ (X,P)X. ^ [- g.cosg-g.sin^.. -j ^^^^. ^^^^^ 

In evaluating the branch line 2 integrals it will be necessary to make some approximations. 
For the frequency range 1 to 1,000 c/s it is permissible to state 

kl ~ coVo^o 

/:!- — ^coMoo•3• (2.5) 

Designating kz=a{\—i), defining a new parameter of integration, t, by 

\=^a{l—i—it) 

and replacing the Hankel function by the first term of its asymptotic series, the branch line 2 
integral for the Hertz potential in region 1 becomes 

I^j,%^-2Ma'^'e'''^%dsJ^ g-«p(i+o (^ t'^2p(t)e-'^klt (2.6) 

\ Trp Jo 

where 

^^~ D{\{t)) 

in which the jS's are expressed in terms of t, and 

The definitions 5i = (7i/a;eo, b^^azlo^ea, ^ = sin fijij and c=cos ^2^ are used. 

The first term of an asymptotic expansion, obtained by integrating terms resulting from a 
Maclaurin series expansion of F{t), gives 

where hn= {kn +2ia^y'^ , n=l, 2, and i==sin 62!!, c = cos 62/^. Since \ki^\>\2ia^\ it is quite accu- 
rate to take hi^h. 

Succeeding terms in the series would go as p"^, p~^, etc. The factor exp [— ap(l+i)], 
or exp [—ikzp], indicates the physical process to be one of propagation along and within the 
ionosphere. 

Similarly the branch line 2 contribution to the Hertz potential in air is found to be 

M2^'^abie~"^''''^e~'^^^^^ ^-apci+o 
Ibl2- !,h fhi:\^h'-o\2 (^1 si^ b2Z-dib2 cos bzz) — ^ (2.8) 

Evaluation of the mode solutions requires knowledge of the poles, X^. Taking reasonable 
values for the physical constants involved, it can be shown that the pole for the lowest order 
mode is given approximately by 
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Approximate formulas for poles of arbitrary order are available (see, for example, Wait 
[1960a and b]). For the frequency range under consideration here only one pole, that which 
for perfectly conducting boundaries would be related to the TEM mode of transmission, is 
very close to the real axis of X. The others are located near and to the right of the negative 
imaginary axis, and their relatively great imaginary part insures that they are heavily damped 
at appreciable distances from the source. 

Thus the important mode solution for the 1 to 1,000 c/s frequency range is, in region 1, 



and in region 2 



27VM Wi^^Hz-zo) 
jj,,__ZwM ^^^^(2)(x,^) ' (2.10) 

2.2. Horizontal Electric Dipole 

The development in this case is similar except that two components, ttx and tt^, of the Hertz 
vector are required to match the boundary conditions, for a dipole oriented along the a:-axis. 
Lack of symmetry in azimuth gives rise to Bessel functions of order 1 in the formulation of the 
TTz component. It is found that for the ic-component 



»xl 



7^x2=^^ 



where Wi(X)= 






=M f" [e-'^>i^+^»i+TFi(X)e'^'<'-^»'] ^ J„(Xp)(/X (2.12) 

Jo 'Pi 

' f" [V2i\)e-'»'''+W2i\)e'^^'] ~ 6-'^''-» J,iXp)d\ (2.13) 

J '^ Pi 

=Mii Vi{X)e-'^''^e-'^'''>Jo{Xp)d\ (2.14) 

J 'iPl 



with 1/2=^2 



|8l+f/2 

/S3 + ii82 tan 02h 



^2+ifi3 tan ^ik 



and V2i\) = -idi ^ 



;S, P2+U2 



iSsft+t^/ 



W2(X) = -i8,i^ 



1 02-U2 



P3P'+U2 



Region 3 wUl not be further considered. 

For the ^-component we have in regions 1 and 2 



where 



,r,i=M cos <t> r 0i(A)e'^i<^-V i^lM x^\ (2.15) 

Jo '2'Pl 

7r,2=M cos <^ i^2(X)e-*^2'+<^2(X)e*V 'j, (Xp) ^-^ \'d\ (2.16) 

Jo L J I iPi L 



. (y, Jl3A{i-&>)il3,s-dil32c)(^2C+iM-5i(i+d3)fi2' 

''^'- ^ {^2C + iM(hfi2C-M(fil+U2)iPl + X2) 
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in which X2 — —i(j2^'2, and 6? and c are defined as before. Contour integration as in the case of 
the vertical electric dipole gives, for the important contributions in region 1, 



Ibl2^''' ^2.^2Mae-''"'G{0)e -^^^ (2.17) 

p 



where 6^(0)- ^ ^ 



with bny n=\, 2, defined as before. In addition 



4^2 _ap(l+i) 

IsL2'"''^Mcos4>-r^,e'H'-''') -2 (2.18) 



and the mode contribution (excited only by ttJ) is 



where 

In region 2 
where 



I^" «2;rM cos 0Xoiy{^^(Xop)e'''"'(^-^») . 1 x ^^^^°^ ' (2.19) 

(__-)(^,+X.)U, 






I'jfll ^2Me""a"'SiO) , — (2.20) 

p 

„ , , _ 2a'^^e~''''"°(6i cos biZ+ibj sin 622) 
^^~ (bic+ib,^' 

I^a^^UI cos 0a^a.e«'/^e-.^o ihcosh.z+ib^^.inb.z) e^^ ^^.21) 

and for the mode we have 

r^^^ =2iTM COS ^\Hi'^ (Xop) ^^^^ .^ ^^V. (2.22) 

where 

^.. .N_ [(^-^i)(fe-^3M(feg+^'M-gi(^+^3)^2](ft sin fe^-gift COS fe^) 
r(\,^)- —-— . 

3. Field Components 
3.1. Vertical Electric Dipole 

The field components are found in the usual way by differentiation of the Hertz potentials. 
In the following the superscript {v,m) will be used to denote 'Vertical electric dipole, mode 
solution^^ (assuming that only the one mode contributes substantially), and (v,b) will denote 
^Vertical electric dipole, branch line solution/^ It is noted that the term (kj—^i^^^) is equal to 
Xq, but otherwise the approximation ^l^^ ^ki can be made throughout. The field components 
in the sea water are : 

E{l-^=-^^^ Xgj?f (Xop) , ' , 7 V > (3.2) 

Hi^^ ™' =2tTMX§^f' (Xop) .jTyy ■> (3.3) 

l.*^ IT" ■^2; >,=>(, 
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Figure 3. Relative magnitude of Eip as it varies with frequency o 

at 1,000 /cm distance from source. Attenuation due to 5 

antenna submergence is not included. uj 

Behavior of other field components is similar, except that the rate of increase of < 

E\z (mode) is greater than that of E\p (mode), as frequency increases. -j 
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In the atmosphere the fiehls are 



-apd + i) 



where 



g(Xo,^) 



TTO), m) . 
-^22 



0-1 
2M7r 



(-l/Xo)(i^i+Z,)Uxo 



0-1 



Xowr(Xop)^-^"^^"^« 



??(Ao,2) 



iy£'-^=-^Xoffr(M^" 



'2* 



5i 



(-l/Xo)(K,+Z,)Ux„ 
P(Xo,g) 



^2°/' - La2-2 ^T2— (^ «os 6^2+5x62 sin 6,0) 

(61 sin 622—5162 cos 622) 



^&« — 



SiSjftp' P- 



5,636^8^ p^ 



, ^4M|aVj^^ ^_,,,,, e^^ ^.^^ 622-5162 cos 622) 
0)^0010302^ P 

q{X,z) = diP2 sin ftg — ft cos fts 

|)(X,2) =5i cos ft2 +|8,/ft sin fts. 
69 



(3.4) 
(3.5) 
(3.6) 



(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 



3.2. Horizontal Electric Dipole 

Using the superscript Qt^m) for the mode solution in this case, and (A, 6) for the branch 
line contribution, it is found that in the sea water, 

£«-«>=2M*i^ COS <^X§[Z?KXop)-Hi(\op)]e«'e-^«) , .nJI^^YV ^^.IS) 

^g.») =4Mt cos ^i^P (Xop)X§e**'^-^o' . ,Awi^iY^> ^3.14) 



H{l-'-'=4TMa, COS <^X^[Fr (Xop)-Sr(Xop)]e"i<^-^o' ._,/. w«^v^^ , (3-15) 



r(Xo) 

(-l/Xo)(^,+Z2)'x=X„ 



where 



Elf''^=cos,p(kj+2ia^)irif + (l-i)kiair','i , (3.16) 

S,<J'" = -sin#Mr, (3.17) 

£'g.6)=cos0^la(l-^)7ri?^ (3.18) 

i7«.«=iLsin<^(M, (3.19) 

COMo 

Z?j<^' " =i^i COS </.ir<?' + (1 +^)aTi?> (3.20) 

H[h. w =sin 0(1 +i)mlf (3.21) 



irff «-2V2aMe«'/4 



^zfci(z-2o) g-ap(l + i) 



(bic+ib2sy p^ 



, ,. 4aWcos0 ,_,^, 62^^)^ 
In the atmosphere the field components are: 

£i^-^=2^7^il/ cos </>X?^?> (M ^ t/tI^ (3-23) 

P2 (P1+A2; X=Xo 



iJ(J-) = _coeoxM cos 0X§[//f (Xop)-^f (Xop)] ^ (^^^^^ - (3.24) 



x=x^ 



where 



.. X [(^-gl)te^-^3M(/32C + ^'fe)-^l(^ + ^3)ffi](igl cos fe^ + 5i/32 sin ig23) 
^ ' ^ (-l/Xo)(/32C + ^/33^)(53ftC-^3^)(^l+C72) 

_[(7;-50(ft^-g3ig2c)(feo+^M-gi(^+^3)ig2](igi sin l32Z-8i^2 cos /322) 



r(X,^) = 



and 



(- 1/Xo) (^2C+iM (83P2C-M (^l+U2) 



£'ip'^^=cos05i(6i cos h2Z+ih2 sin 62^) ^ — (3.25) 

P 

p-ap{l+i) 

£'iS'^^ =— sin (pAiklibi cos 62^+^62 sin M ^ — (3.26) 

p 

^-ap(l+z) 

£'^J. 6)= COS 052(62 COS b22+ibi sin M) 2 — (3-27) 

P 
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i7g*^^ = -sin <A ^^^^^' (61 sin b2Z-ib2 cos 62^) ^- (3.28) 

H^4'''^=Gos<j>B^{b2 COS 62^+^61 sinM ^ — 2— (3-29) 



-apa+i) 

^?2^J•*^=sin0^la(l + ^)(6l cos 62^+^62 sin M ^—, (3.30) 

P 



i:^^^'•"^=sln0^la^-t-^j(^Ol cos (?2^i-'?'^2 sin 02^; 
where 



^1= 






B,= {kl+2id')A,+a{l+i)b2A2, 

B2=klA2+a(l-i)b2A,, 

53=a(l+^)^2+^Ml. 

4. Conclusions 

Not too many numerical results have been obtained from the formulas as yet. Some calcu- 
lations have been carried out, however, for the fields in sea water at a distance of 1,000 km and 
frequencies of 1, 10, 100, and 1,000 c/s. Conductivities of 4 and 10"^ mhos/m were assumed for 
sea water and ionosphere, respectively. The ionosphere height was taken to be 90 km, a figure 
which has been used with some success in ELF studies [Wait and Carter, I960]. 

The more important results can be summarized as follows: (See also fig. 3.) 

1. The mode solution (corresponding approximately to the TEM mode in a parallel-plane 
waveguide) involved here is excited more efficiently by the vertical than by the horizontal 
electric dipole. This is not unexpected. 

2. The branch line contribution (corresponding to energy transmitted along the air- 
ionosphere boundary, but within the ionosphere) in the case of the horizontal dipole is greater 
than the mode contribution at 1 and 10 c/s by as much as 10^ depending on the component 
concerned. 

3. At 1 c/s the field components due to the horizontal dipole are greater than those of a 
vertical dipole of the same moment by virtue of the branch line contribution. At 10 c/s and 
higher the field components of the vertical dipole are greater by virtue of the mode solution. 
At greater ranges, the vertical dipole mode solution would predominate even at 1 c/s, due to 
its much smaller attenuation with distance. 

4. The system attenuation between half wavelength (in water) ^^ coaxial antennas'' [Moore, 
1951] separated by 1,000 km is more than 250 decibels for any frequency in the range consid- 
ered, not including attenuation in the ^-direction, that is, attenuation due strictly to the sea 
water. In view of the fact that noise power is substantial in this frequency range, the difficulty 
of practical communication under such conditions would seem to be great. 

5. At distances less than 1,000 km the branch line contribution will be important even in 
the higher portion of the 1 to 1,000 c/s frequency range, and should therefore be investigated 
more carefully. Since the asymptotic series expansion is usable only for large values of p, 
some method other than that of asymptotic expansion must be used. Perhaps numerical 
integration would be in order. At closer distances it is possible that substantial contributions 
from the more highly damped modes may occur; therefore these may also require investigation. 

More details as to the actual figures obtained by computation are available in a report, by 
the writer, published by the University of New Mexico, Engineering Experiment Station, 
Technical Report EE-44, February, 1961. 
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The reader will note by reference to figure 3 that since the magnitude of the field compo- 
nents increases with frequency in the case of the mode solutions, and since attenuation due to 
travel in sea water increases monotonically with frequency, then at any combined depth of 
receiving and transmitting antennas there is evidently an optimum frequency. There is no 
optimum for the branch line solutions, however. 

It should be mentioned that in view of the large system attenuation at distances of 1,000 
km or greater, it is perhaps only of somewhat academic interest to consider refinements in the 
model which would take into account the anisotropy and inhomogeneity of the ionosphere. 
As pointed out by Wait [1960a] the assumption of a refractive index increasing exponentially 
with height causes an attenuation rate which varies more rapidly with frequency than that for 
a homogeneous ionosphere. The experimental data seems to require this more elaborate 
model. In the problem considered here, however, experimental testing of the theory is prob- 
ably not feasible except at distances less than 1,000 km or so, and at such distances it is not 
expected that the results will be drastically changed by the different attenuation rate. The 
effect of the anisotropy caused by the earth's magnetic field is likewise expected to be relatively 
minor [Wait, 1960a] at such distances. 
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